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ABSTRACT 

Non-thermal pressure from turbulence and bulk flows is a fundamental ingredient in 
hot gaseous halos, and in the intracluster medium it will be measured through emission 
line kinematics with calorimeters on future X-ray spacecraft. In this paper we present 
a complementary method for measuring these effects, using forbidden FUV emission 
lines of highly ionized Iron which trace 10^ K gas. The brightest of these is [Fexxi] 
A1354.1. We search for these lines in archival HST-COS spectra from the well-known 
elliptical galaxies M87 and NGC4696, which harbor large reservoirs of 10^ K gas. 
We report a 2.2a feature which we attribute to [Fexxi] from a filament in M87, and 
positive residuals in the nuclei of M87 and NGC4696, for which the 90% upper limits 
on the line flux are close to the predicted fluxes based on X-ray observations. In a 
newer reduction of the data from the Hubble Spectroscopic Legacy Archive, these 
limits become tighter and the [Fexxi] feature reaches a formal significance of 5.3cr, 
neglecting uncertainty in fitting the continuum. Using our constraints, we perform 
emission measure analysis, constraining the characteristic path length and column 
density of the ^ 10^ K gas. We also examine several sightlines towards filaments 
or cooling flows in other galaxy clusters, for which the fraction of gas at 10^ K is 
unknown, and place upper limits on its emission measure in each case. A medium- 
resolution HST-COS observation of the M87 filament for ~10 orbits would confirm 
our detection of [Fexxi] and measure its width. 

Key words: galaxies: clusters: intracluster medium; galaxies: haloes; galaxies: indi¬ 
vidual: M87; galaxies: individual: NGC 4696; ultraviolet: galaxies 


1 INTRODUCTION 

Turbulence and bulk flows provide non-thermal pressure in 
the hot gaseous halos of galaxies, galaxy groups, and galaxy 
clusters. On large scales, this is predicted by numerical sim¬ 
ulations (e.g. [Norman &; Bryan||199^|Frenk et al.|1999| ) due 
to mergers and the accretion of satellites. In the centers of 
these systems, turbulence and bulk flows are also introduced 
by active galactic nucleus (AGN) feedback (jChurazov et al.j 


20011; other mechanisms such as cosmic ray-driven instabil 
ities (Sharma et al. 20091 and galactic stirring ( Kim|[2007 


jRuszkowski fc Oh||2010| ) may also be important in the intr¬ 
acluster medium. 

These motions play a fundamental role in the physics 
of these systems. They transport energy, mass, and met- 
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als from the central nucleus into the diffuse halo, affecting 
the thermodynamic state of the intracluster medium and 
mediating the enrichment of the intergalactic medium. The 
non-thermal pressure from turbulence is also necessary to 
understand in order to estimate the masses of galaxy groups 
and clusters from hydrostatic equilibrium arguments. Large- 
scale flows are expected as well around the supermassive 
black holes at the centers of these systems, particularly as 
one approaches the Bondi radius | Bondi|[i952 K In order to 


understand turbulence and bulk flows in better detail, 
great deal of effort has been focused on the kinematics of 
the hot halo gas in galaxies, galaxy groups, and galaxy clus- 


et al.|2004l|Schuecker et al.|2004l|Werner et al.|2009| Sanders 

et al.|2011l, and simulations (e.g. Evrard|1990 

Schindler & 

Mueller|1993| Norman & Bryan|1999| Nagai et a 

l.|2003 

|Sun- 

yaev et al.||20031 Briiggen et al.||2005 [Briiggen et af. 
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2 Anderson and Sunyaev 


'Va.zza, et al.||200^ |Lau et al.||2009[ [Battaglia et al.||2012[ 
Gaspari et al. |2014[ [Zhuravleva et al.||2014|). 

However, these measurements are extremely challenging 
and expensive. Constraints can be obtained in some cases 
through the kinetic SZ effect ( Sunyaev fc Zeldovich|[T980 1, 
and observations of resonant scattering in the intracluster 
medium (Gilfanov et al. 19871. The more direct approach 


is to measure the broadening and redshift of emission lines, 
which trace turbulence and bulk motions respectively. For 


bulk flows, Inogamov & Sunyaev (20031 showed that the 


typical motions of bulk flows in galaxy clusters are much 
larger than the peculiar velocities of the clusters, making 
such measurements possible. Also, the thermal velocities of 
heavy ions in hot gas (uth = \/2kT /Mi) are ~ yw rup 
times smaller than the speed of sound (cs = 'ykT/fimp) 
in the same gas. Measuring the broadening of (for example) 
Iron lines in hot plasmas, we can detect broadening due to 
turbulence because even subsonic turbulence with velocities 
on the level of 10% of sound speed increases significantly 


the width of heavy ion lines (Inogamov & Sunyaev 20031. 


To do this, one needs to resolve the emission lines, which are 
primarily the soft X-ray Ovill and Fe L-shell lines for gas 
around 10^ K, but CCD detectors on modern X-ray tele¬ 
scopes (R = A/AA ~ 15 — 20) lack the required spectral 
resolution. X-ray grating spectrographs have better resolu¬ 
tion {R ~ 100 — 1000) and can identify important features 
like the Fe L shell lines, but these lines lie very close to 
each other and are usually blended. Using a calorimeter, the 
Hitomi satellite was recently able to make a high-precision 
measurement of the broadening in the Fe Ka line for a sin¬ 
gle target - the Perseus cluster (Fabian et al 2016), which 
is the X-ray brightest cluster in the sky. In coming decades, 
such devices aboard next-generation spacecraft like Athena 
or the X-ray Surveyor will offer similar improvements in clus¬ 
ters hot enough to produce this line. For 10^ K gas, how¬ 
ever, X-ray calorimeters will be limited to R~300-400 for 
the prominent X-ray emission lines. 

Here we discuss a complementary approach: measuring 
the kinematics of far-ultraviolet emission lines. In the FUV, 


the Gosmic Origins Spectrograph (COS) instrument (Green 


et al. 20121 on the Hubble Space Telescope (HST), offers 


much higher throughput and spectral resolution. Typically 
ultraviolet lines are characteristic of cooler plasma, at tem¬ 
peratures of 10“^ — 10® K, but there are exceptions. Some 
FUV lines have been used as diagnostics of plasma at tem¬ 
peratures around 10^ K, both in Solar flare spectra (e.g. 
Doschek^^etalJ[19^ Chenge^^ay[1979|^ Masonet^^alJ[1986) 


Feldman & Doschek 1991[ Innes et al. 2003| Innes et ah] 

2003 Doschek & Feldman 2010[ Young et al. 20151, and 


in the coronae of nearby late-type stars (e.g. Maran et al. 


1994 


2003 


[Linsky et al.||1998[ [Ayres et al.||2003j [Redfield et al. 


Dupree et al. 12005 I. In this paper we will discuss these 


lines in an extragalactic context, as tracers of the intraclus¬ 
ter medium and accretion flows. COS has one low-resolution 
grating (G140L) with R ~ 10® and two medium-resolution 
gratings (G130M and G160M) with 7? ~ 1 —2 x 10^, all three 
of which we utilize in this work. We note that this method 
is conceptually similar to a method proposed by jSyunyaevj 
& Churazov (1984 \ to use a 3 mm hyperfine transition from 
Lithium-like ®'Fe. 

In Figure 1 we present basic diagnostic properties for six 
of the brightest of these lines. They are all ground-state mag¬ 


netic dipole transitions of various species of highly ionized 
Iron. In the top panel, we compare the fractional abundance 
of these six species of Iron for plasma in collisional ioniza¬ 
tion equilibrium (GIE), based on the CHIANTI database 
V. 7.15 (Dere et al. |1997| Landi et al. 20131. In the mid¬ 
dle panel, we present the effective cross section for absorp¬ 
tion from these lines, scaled to the Thomson cross section 
(6.65 X 10“®® cm“®). Only thermal broadening is considered, 
and we assu me an Iron abundance r elative to Hydrogen of 
3.16 X 10“® ( Grevesse fc Sauval|l998 1. In most extragalactic 
cases these lines will be optically thin, and since they are 
forbidden lines the optical depths will generally be very low. 
They may be detectable in some situations, although in this 
paper we focus on the lines in emission. It is important to 
emphasize that, since these lines are forbidden transitions 
and optically thin, the observed line widths are not affected 
by most radiative transfer effects which could otherwise in¬ 
troduce artificial broadening to emission lines. 

In the bottom panel of Figure 1 we present the emissiv- 
ity of these lines, defined such that e = LI{Zpi,J nf.npdV) 
where L is the line luminosity and Zpe is the fractional 
abundance of Iron relative to Hydrogen, for which we as¬ 


sume the Solar value measured by Grevesse & Sauval (19981. 


The quantity f neripdV is the volume emission measure, 
which for a pencil-beam sightline differs from the pencil- 
beam emission measure f Uf-ripdl by a factor of the area A 
covered by the sightline. This emissivity is also derived from 
the CHIANTI database, by multiplying their G(T) function 
from gofnt. pro by dvr steradians. It is computed assum¬ 
ing collisional ionization equilibrium and that the lines are 
optically thin. 

In Figure 1, the brightest of these lines is [Fexxi] 
A1354.08. In this paper we will focus primarily on this line, 
though in some cases information on this line is not available 
and we consider other lines. These lines trace plasma in a 
narrow range around 10^ K, with the emissivity of [Fexxi] 
peaking at 10^ ®® K. The instrument of choice for emission¬ 
line searches is COS, since it has the highest throughput 
in the FUV at present. We searched the COS archive for 
observations of systems known to host a significant reser¬ 
voir of gas around 10^ K. We considered M87, Cen A, NGC 
4696, and HCG 92; out of these targets, M87 and NGC 4696 
have already been observed by COS in the FUV. There are 
also observations of HCG 92, but the the sightlines target 
a shock front near NGC 7318b which is likely significantly 
cooler than 10^ K. We also searched for observations of Ha 
filaments and/or cooling flows in galaxy clusters; the ambi¬ 
ent media here are hotter than 10^ K, but at the locations of 
these features multiphase gas is present, which may include 
10^ K gas. We find suitable observations for A1795, Perseus, 
Phoenix (SPT-CL J2344-4243), and Zw 3146. In these latter 
cases we lack definite predictions of massive reservoirs of 10^ 
K gas along these sightlines, but we can use upper limits to 
establish independent constraints on the differential emis¬ 
sion measure. A summary of the observations we consider is 
given in Table 1. 

In the next several sections of this paper, we describe 
the results of our search. In Section 2, we present a ten¬ 
tative detection of [Fexxi] and weak evidence for [Fexix] 
in the filaments projected 1.9 kpc from the center of M87. 
In sections 3 and 4 we present upper limits on [Fexxi] in 
the centers of M87 and NGC 4696. In Section 5 we dis- 
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Figure 1. (top) Fractional abundance of six species of highly ionized Iron, for plasma in CIE. (middle) Effective absorption cross sections 
for the strongest FUV transition for each species, scaled to the Thomson cross section and assuming the lines are optically thin and 
affected only by thermal broadening, (bottom) Approximate emissivity for these transitions (e = L/(ZFe f nenpdV)), in units of erg cm^ 
s~^. Middle and bottom panels assume Solar abundance. 
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cuss potential sources of contamination, i.e. other emission 
lines which could be confused with [Fexxi] and/or [Fexix] 
in these galaxies. In Section 6, we report limits on [Fexvill], 
[Fexix], [Fexxi], and [Fexxiii] in COS observations of sev¬ 
eral additional clusters of galaxies. In Section 7, we show 
how the results change if the co-addition and error estima¬ 
tion is modified. Finally, in Section 8 we discuss prospects 
for future searches. In the Appendix we discuss the filtering 
for geocoronal emission performed for the G140L spectra 
towards M87. 


1.1 Observing Extended Emission with COS 

It is important to emphasize that COS is not optimized for 
the analysis of extended emission. It employs a circular slit¬ 
less aperture, so if emission is extended in the sky along the 
COS dispersion axis, the lines will appear artificially broad¬ 
ened in the resulting spectrum ( Green et al.|2012 l. This ef¬ 
fect is likely contributing to the abnormally large line widths 
seen for Civ and Hen, and is also plausibly occurring at the 
locations of [Fexxi] and [Fexix]. 

To estimate this effect, we fit a Gaussian to the along- 
dispersion vignetting profile given in Table 6 of the COS 
instrument science report 2013-03 ( Proffitt et al. |2010| and 
use the along-dispersion plate scale of 0.0285”/pixel to con¬ 
vert into velocity units at 1354.1A. From this analysis, we 
estimate a — 648 km/s for G140L and cr = 81 km/s for 
G160M. We consider these values the theoretical prediction 
for spatial broadening of uniformly distributed [Fe xxi] emis¬ 
sion. They are a bit lower than the a ~ 110 km/s estimated 
for spatially broadened lines using G130M observations of 
the compact blue galaxy KISSR242 ( [France et al.|[2010[ ), 
but the latter measurements include the intrinsic line widths 
from the galaxy as well. 

An additional effect is that calcos pipeline is intended 
for the analysis of point sources, so the fluxes inferred for ex¬ 
tended emission like [Fe xxi] may be biased. To estimate this 
effect, we fit a fifth-order polynomial to the cross-dispersion 
profile from Proffitt et al. I polot (their Table 5) and inte¬ 
grate this profile over the 2.”5 aperture in order to estimate 
the expected throughput for a well-centered point source. 
We find that the expected throughput for uniformly dis¬ 
tributed extended emission is about 85% of the expected 
throughput for a point source. However, since the exact cor¬ 
rection depends on the (unknown) spatial distribution of the 
emission, and the correction is much smaller than the cur¬ 
rent measurement uncertainties, we neglect it in the present 
analysis. 


to be bright in C iv A1549 (Sparks et al. 20091. A subsequent 


observation with the COS G140L grating shows unambigu- 


2 M87 FILAMENT 

M87 is the central dominant elliptical galaxy at the core of 
the Virgo cluster, and is suffused with hot X-ray emitting 
gas. Near the center of the galaxy, much of this hot gas is 


multiphase (Belsole et al. 2001 Molendi 2002 Simionescu 


et al.|2008 Werner et al.|20l0 l, spanning at least a factor of 


six in temperature from about 3 keV to about 0.5 keV, with 
absorption possibly obscuring gas below this temperature 


ous e mission lines from both C IV and the He ll (Sparks et al. 
20121, indicating the presence of 10® K gas in these struc¬ 


tures. At lower temperatures, Werner et al.| ( [2bl3 l used the 
PACS instrument on Herschel to detect [Cii] A158/im emis¬ 
sion from M87 which seems to trace the filamentary struc¬ 
ture as well; this is typically the dominant cooling line when 
Hydrogen is neutral ( Dalgarno fc McCray|1972 \ and is most 
effective at temperatures around lO'^ K. 

Low-resolution G140L observations of one of the FUV- 
bright filaments are available, which we examine here. We 
downloaded and examined the xldsum files for datasets 
LBK104010, LC7T01010, and LC7T01010, which contain 
the merged spectrum as reduced by the calcos pipeline. At 
the redshift of M87, the geocoronal O l A1356 falls very close 
to the expected location of [Fexxi] A1354.1, so we first con¬ 
struct a spectrum from the data taken only from orbital 
night (using 7858 s out of a total exposure time of 31798 
s) in order to exclude the possibility of geocoronal emission 
contaminating our search for [Fexxi]. We discuss this pro¬ 
cess in more detail in the Appendix. 

We generate a series of bins approximately 14 pixels in 
size, and within each bin we average together the fluxes from 
each pixel of each spectrum that falls within the bin. A por¬ 
tion of the night-only spectrum from segment A is shown in 
Figure 2. The FWHM of the line spread function is about 
7 pixels, so this binning combines about two resolution ele¬ 
ments (resels) per bin. We use the calcos pipeline errors for 
the error spectrum. We use a redshift of u = 1284 ± 5 km/s 
( Cappellari et al.|2011 l to estimate rest-frame wavelengths. 

The Hell multiplet and the Civ doublet, which trace 
10® K gas, are clearly visible and have already been dis¬ 
cussed by Sparks et al. ( 2012| ). A weak positive residual is 
also visible at the expected location of [Fexxi]. We perform 
Gaussian fits to the observed [Fexxi] line profile, after sub¬ 
tracting a linear continuum fit to regions of the spectrum 
outside the line (Figure 2, lower left panel). We estimate 
likelihoods over a grid of three parameters (mean velocity 
Wo, velocity dispersion a, and integrated flux F) using the 
statistic, and compute medians and central 68% confi¬ 
dence intervals for each parameter, marginalizing over the 
other two. The best-fit profile has wo = —250 km/s, u = 356 
km/s, and F = 7.2 x 10“^^ erg s“^ cm“^. The best-fit profile 
is shown in Figure 2 (lower right panel). However, the sta¬ 
tistical significance of this detection is only 1.0a, so it is not 
formally significant. The 90% upper limit on the integrated 
flux is 1.7 X 10“^® erg s“^ cm“^. Assuming a distance to 



M87 of 16.7 Mpc (Blakeslee et al. 20091, this corresponds 
to an [Fexxi] A1354.1 line luminosity of L < 5.8 x 10®® erg 
s“® within our pencil-beam aperture (or if the emission is 
uniformly filling the field, this limit increases by about 15%; 
see section 1.1). 

We next tested the strength of the [Fexxi] feature as 
a function of the solar altitude relative to the instrument. 
We found that, with the exception of three observations in 
dataset LC7T02010 (which are discussed in the Appendix), 
the line flux remains constant, implying that Ol A1356 air- 
glow emission is negligible in these spectra (this is also true 
for the M87 nuclear spectra discussed in the next section). 
We therefore also consider the combination of all the M87 
filament data (daytime and night-time, excepting the three 
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Table 1. Observations Examined in this Work 


Dataset 

target 

Grating 

Central Wavelength 

(A) 

Exposure Time 

(s) 

LBK104010 

M87 filament 

G140L 

1280 

5025 

LC7T01010 

M87 filament 

G140L 

1280 

13386 

LC7T02010 

M87 filament 

G140L 

1280 

13387 

LC8701010 

M87 nucleus 

G130M 

1318 

1000 

LC8701020 

M87 nucleus 

G130M 

1309 

1000 

LC8701030 

M87 nucleus 

G130M 

1300 

1200 

LC8701040 

M87 nucleus 

G130M 

1291 

1200 

LC8702010 

NGC 4696 nucleus 

G130M 

1291 

1191 

LC8702020 

NGG 4696 nucleus 

G130M 

1300 

1191 

LC8702030 

NGG 4696 nucleus 

G130M 

1318 

1417 

LC8702040 

NGG 4696 nucleus 

G130M 

1327 

1417 

LClJOlOlO 

A1795 filament 

G140L 

1280 

16443 

LBPOOlOlO 

Perseus filament 

G140L 

1105 

11437 

LBP002010 

Perseus filament 

G140L 

1280 

8310 

LC9552010 

Phoenix nucleus 

G160M 

1577 

19752 

LCDI03010 

Zw 3146 nucleus 

G140L 

1280 

25279 

LCAY02010 - LCAY06040 

HCG 92 

G130M / G160M 

several 

42074 


List of the HST-COS observations examined in this work, along with basic properties of 
the observations. For HCG 92, we have combined together 20 observations of five different 
regions along the same shock front. 


contaminated observations mentioned above). This spec¬ 
trum is shown in Figure 3. For this spectrum, we measure 
for [Fexxi] best-fit line parameters of vo = —3901^4° km/s, 
a — 4001i 20 km/s, and F — T.Slg g x 10“^^ erg s“^ cm“^. 
These parameters are nearly identical to the night-only val¬ 
ues. For Ol A1356 with an expected 3:1 ratio between the 
lines in the doublet, the velocity relative to the position of 
[Fexxi] is about —820 km/s, so this line is about halfway 
between the expected locations for [Fe xxi] and airglow emis¬ 
sion. However, the night-only line also has a similar negative 
velocity, and there are additional uncertainties such as a 150 
km/s uncertainty in absolute wavelength calibration, so we 
believe the line can properly be attributed to [Fexxi]. The 
negative velocity is a bit of a concern, however, and prob¬ 
ably unphysical; follow-up observations at better spectral 
resolution (see section 8) will hopefully help clarify the sit¬ 
uation. As mentioned above, the flux is nearly identical to 
the night-only case, but the longer integration time yields a 
larger S/N, so that the significance of the feature increases 
to 2.2cr in this spectrum. 

This spectrum also contains information about [Fe xix] 
A1118.1. In Figure 4 we display the merged G140L spectrum 
of segment B, which covers the shorter-wavelength portion 
of the FUV, using all the observations from all three datasets 
targeting the filament region, as in Figure 3. Below llOOA 
the sensitivity of COS declines dramatically, and at 1134A 
NI airglow emission is visible, but between these two re¬ 
gions the [Fexix] A1118.1 line can be constrained. There is 
a weak feature at this location, significant at just 0.3 ct, with 
the best-fitting Gaussian having vo = —190 km/s, a = 346 
km/s, and F = 1.6 x 10“^^ erg s“^ cm“^. The 90% upper 
limit on the integrated flux is 8.6 x 10“^^ erg s“^ cm“^, 
corresponding to a line luminosity L < 2.9 x 10^® erg s“^. 

In Figure 5, we illustrate an emission measure analy¬ 
sis based on these constraints. Using the emissivities shown 
in Figure 1, we derive the pencil-beam emission measure 
f UeUpdl = L/{ex Ax Zpe) corresponding to the constraints 


on the line luminosities, where A is the physical area sub¬ 
tended by our aperture at the distance of M87 (3.0x 10^ pc^) 
and Zpe is the fractional abundance of Iron relative to Hy¬ 
drogen in the hot plasma, which we assume to be 3.16 x 10“® 
(i.e. Solar abundance; [Grevesse fc Sauval|199^ . We consider 
the 90% upper limit for [Fexix] and the 2.2a measurement 
for [Fexxi], although we also show the 90% upper limit on 
[Fe xxi] from the night-only spectrum in order to be conser¬ 
vative. 


The resulting curves for [Fexxi] have a clear minimum 
around T = 10^'®® K, which is where the emissivity of 
[Fexxi] is highest. This is the portion of the plot at which 
one would naturally expect physical solutions.The implied 
emission measure at this temperature is 1.4±0.6x 10^® cm“®. 
[Fexix] has a lower emissivity than [Fexxi], so it does not 
give as strong of a constraint; the 90% upper limit on the 
emission measure at its peak temperature is 5.9 x 10^® cm“®. 


Next we estimate the physical properties of the hot gas 
in our pencil-beam aperture. To do this, we combine our con¬ 
straints with an additional piece of information: the pressure 
of the gas as measured from X-ray observations. We adopt 
a value for the pressure of P = 1.91Pe = 4 x 10® cm“® K, 
based on Churazov et al. (20081. This value is fairly robust: 


Di Matteo et al. (20031, Ghurazov et al. (20081, and Rus¬ 


sell et al. (20151 all find similar results, and they all find 


that the deprojected pressure profile is approximately con¬ 
stant within the central 2.5 kpc, except for an increase in 
the central 100 pc. 


Combining our results for the emission measure with 
the X-ray measurements of the pressure, we can derive sev¬ 
eral additional constraints. The equivalent electron density 
is rig = PelT, which we plot in the upper axis of Figure 5. 
The corresponding electron column in our aperture is given 
by 
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Figure 2. (top) Average spectrum of COS G140L observations of a filament projected 1.9 kpc from the nucleus of M87. The spectrum 
is shown in black and the Icr error per bin is shown in green. The spectrum has been binned to a resolution of 14 pixels (approximately 
two resels). Prominent airglow lines are indicated in light blue. Since Ol A1356 falls very close to the expected location of [Fexxi], we 
have filtered this observation to display only the spectrum taken during orbital night, which should eliminate any detectable emission 
from this airglow line (see Appendix). There is a weak suggestion of [Fexxi] in this spectrum, but it is not formally significant due to 
the large errors per bin. The lower left panel shows the spectrum within 2500 km/s of [Fexxi], as well as our linear fit to the continuum. 
The lower right panel shows the background-subtracted spectrum within 2000 km/s of [Fexxi] as well as the best-fit Gaussian profile, 
which corresponds to a l.Oir detection. The 90% upper limit on the integrated [Fexxi] A1354.1 flux is 1.7 X 10“^® erg s"'^ cm“^. 


Ne = Uel = 1.6 X cm ^ X 


EM 


2.5 X 1020 cm -5 


107.05 K 


4 X 10® cm ® K 




( 1 ) 


where the coefficient is derived by assuming p — 1.91pe and 
a proton-to-electron ratio of 0.83. Note that this implies that 
[Fexxi] and [Fexix] would indeed be optically thin, as we 
had assumed when deriving the emissivity. The correspond¬ 
ing effective path length of [Fexxi]-emitting gas in our aper¬ 
ture is 


I = 3.0 kpc X 


X 


T 


107.05 K 


EM 


2.5 X 1020 cm“® 


4 X 10® cm”® K 



( 2 ) 


M = pempNsA, where A is the physical area subtended 
by our aperture at the distance of M87 and we assume 
Pe = 1.2. The metallicity dependence enters through the 
emission measure, which scales as {Z/Zq)~^; if we assume 
the pressure is robustly determined, then the path length, 
column density, and hot gas mass all scale as {Z/Zq)~^, as 
shown in equations and 

At 10^'®® K, our estimates for these quantities are 
Ns = 9 ± 4 X 10^® cm“^, I = 1.5 ± 0.7 kpc, and M = 
2.8 ± 1.2 X 1O®M0. These values are physically plausible for 
the gas around 10^ K in M87, and they also agree with X- 
ray measurements. For example, Di Matteo et al. (20031, 


Churazov et al. (20081, and Russell et al. (20151 all find de- 


projected electron density profiles which are fairly constant 
from 100 pc out to a few kpc, at values of 0.1-0.2 cm~®. Ad¬ 
ditionally, the sightline is projected 1.9 kpc from the center 
of M87, so an effective path length of this order is a reason¬ 
able first-order assumption. 


We use these relations to plot curves indicating Ns and I There is also an X-ray measurement of the emission 

in Figure 5. The total mass in hot gas traced by [Fexxi] is measure of 1 keV gas, which is an independent measure- 
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Figure 3. (top) Average spectrum of COS G140L observations of a filament projected 1.9 kpc from the nucleus of M87. Unlike the 
previous Figure, we have not restricted this spectrum to the night-only observations; O I A1304 airglow is therefore extremely prominent, 
but we argue that the feature near Ol A1356 is probably [Fexxi], with a negligible airglow contribution. The feature is significant at 
2 . 20 -. 


ment against which we can compare our result. This mea¬ 
surement comes from the 1 keV component of the three- 
temperature spectral modeling for this region performed 
by Werner et al. (20131. From their high-resolution Chan¬ 
dra map, they estimate the emission measure of 1 keV gas 
from the sightline corresponding to our observation to be 
around 2 — 4 x 10®° cm“® arcsec”^, using a metallicity 
of 1.46.^0 relative to the abundance tables of|Grevesse fc| 


Sauval (19981. This corresponds to an [Fexxi] A1354.1 flux 
of 3 — 5 X 1O“^^(Z/Z0) erg s“^ cm“^. This is within la of 
our estimate of 7.8 ± 3.6 x 10“^^ erg s“^ cm~^, assuming 
Solar abundance. 


2.1 Civ and He II 

It is interesting to compare the lines tracing 10^ K gas to the 
stronger C iv and He ll lines which are characteristic of 10® 
K gas in this spectrum. We fit a Gaussian to Hell A1640.4 
and a double Gaussian (with normalizations fixed to the 
expected 2:1 ratio) to Civ A1548.2 and Civ A1550.8 in the 
full G140L filament spectrum (i.e. without the filtering to 
exclude data taken during orbital day). We measure velocity 
dispersions of around 600 km/s for both lines, which are 
far larger than the expected thermal velocities for C and 


He ions in 10® K gas (10 and 20 km/s respectively). We 
interpret these dispersions as effects of spatial broadening, 
as discussed in the next section. 


Both lines also appear redshifted relative to the veloc¬ 
ity of M87, by about 150 km/s. The absolute wavelength 
calibration of G140L is 150 km/s, so this apparent redshift 
could be explained by a wavelength miscalibration, but we 
think it is likely to be real. This filament has also been ob¬ 
served in Ha ( [Sparks et al.|199^ and [C ll] A158jrm ( jWernerj 
et al.|[2013 1, and both observations show the same redshift 


for the filament. This bulk velocity corresponds to ~40% of 
the sound speed of the hot 10^ K ambient medium, which 
is evidence that bulk motions through the hot medium can 
be considerable. Interestingly, the weak [Fexxi] feature and 
the weaker [Fe xix] feature both show no sign of redshift rel¬ 
ative to M87; if these lines are confirmed, this would imply 
the 10^ K gas is not moving relative to the galaxy (i.e. the 
filaments are flowing through the hot medium). Deeper ob¬ 
servations with higher spectral resolution are necessary to 
confirm this, however. 


Finally, we consider the fluxes of these lines. We mea¬ 
sure integrated fluxes of 3.6 ±0.9 x 10“^® erg cm“^ s“^ and 
3.6 ± 2.0 X 10“^® erg s“^ cm“^ respectively for Civ and 
Hen, which are roughly consistent with the measurements 
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Figure 4. (top) Average spectrum of the shorter-wavelength segment B for three COS G140L observations of the a filament projected 
1.9 kpc from the nucleus of M87. The spectrum is shown in black and the Icr error per bin is shown in green. The spectrum has been 
binned to a resolution of 14 pixels (approximately two resels). Prominent airglow lines are indicated in light blue. The sensitivity of COS 
declines significantly below llOOA, so the spectrum is entirely noise-dominated below this wavelength. There is a small positive residual 
at the location of [Fexix] in this spectrum, which is only significant at 0.3 (t. The lower left panel shows the spectrum within 2500 km/s 
of [Fexix], as well as our linear fit to the continuum. The lower right panel shows the background-subtracted spectrum within 2000 km/s 
of [Fexix] as well as the best-fit Gaussian profile, which corresponds to a O.Scr detection. The 90% upper limit on the integrated [Fexix] 
A1118.1 flux is 0.9 X 10“^® erg s“^ cm“^. 


of Sparks et al. (20121. In Figure 6 we show the CIE emis- 
sivities for these two lines. They both peak at around 10® K, 
but for Solar abundances their expected ratio at this tem¬ 
perature is about 110, much higher than the observed ratio 
of 1.0. This places a constraint on the shape of the differ¬ 
ential emission measure (DEM) in this region: it must have 
far more weight towards gas at temperatures of 2 x 10® K 
and above, as compared to 10® K. 

Thermal conduction is a likely explanation for this un¬ 
usual line ratio. A cooling flow is an alternative hypothesis, 
but a cooling should have a DEM which is inversely propor¬ 


tional to the cooling function (Fabian 19941, which could 


only reduce the C iv / He ll ratio by an order of magnitude; 
this is insufficient to match the observed ratio. The thermal 


conduction model explored by Sparks et al. (20091 also sig¬ 


nificantly underpredicts Hell, although their updated model 
has more success (Sparks et al. 20121. Additionally, collisions 


with the hot electrons fro m the ambient medium cou ld sup¬ 
press the fraction of (Hicks & Canizares 20011, while 


simultaneously boosting Hell since the latter is composed 
of strong recombination lines (it is the Ha transition for 


Hydrogen-like Helium). Whether or not the latter effect is 
necessary, the physical picture we infer is that C iv traces 
the 10® K gas, while Hen has contributions from a wide 
range of temperatures. 

It is worth mentioning the possibility of photoioniza¬ 
tion being important as well, especially since photoioniza¬ 


tion commonly produces C iv / He ll ratios near unity (Allen 


et al. 19981. While the radiation field from stars and the 


hot gas in M87 is insufficient to affect the ionization bal¬ 
ance significantly here, perhaps a burst of activity from the 
AGN in the recent past might have contributed the neces¬ 
sary photons. This underscores the difficulty in interpret¬ 
ing the filamentary line emission. Clearly the exact physical 
mechanisms producing the C iv and He ll are not yet fully 
understood, complicating a kinematic analysis of this region. 


3 M87 NUCLEUS 

A sightline with the medium-resolution G130M grating tar¬ 
gets the central region of M87, including the active nucleus 


© 0000 RAS, MNRAS 000, 000-000 























































































































9 


[cm 

0.25 0.2 0.15 



T X 10^ K 


Figure 5. Constraints on the pencil-beam emission measure 
f UeUpdl from the filament sightline towards M87 inferred from 
our flux measurements for [Fe xxi] (solid curve and shaded region 
for the measured value and Icr uncertainties for the spectrum 
in Figure 3; the dashed curve shows the 90% upper limit from 
the night-only spectrum for comparison) and [Fexix] (upper left 
curve). We also compute the corresponding average electron den¬ 
sity (upper X-axis), electron column density (dotted lines), and 
path length (dashed lines) for the hot gas, assuming a constant 
total pressure of P = 1.91Pe = 4 x 10® cm ^ K in the sight¬ 
line ( [Churazov et all|2008| |. [Fexix] has a lower emissivity than 
[Fexxi], so it does not yield as strong of a constraint. The [Fexxi] 
curve is sharply peaked at 10^ '^® K, and if the [Fe xxi]-emitting 
gas has this temperature, we obtain physically plausible parame¬ 
ters for the hot gas in the sightline. The shaded dark rectangular 
region indicates the measured emission measure of 1 keV gas from 
the X-ray observations of [Werner et al.| l |2013| l , which we have 
scaled to an assumed Solar metallicity. This measurement is con¬ 
sistent with our estimate of the [Fexxi] luminosity. The physical 
values in this plot are all computed assuming Solar abundance, 
and scale with metallicity as discussed in the text. 


at the edge of the aperture but excluding the jet and the 
bright knot HST-1. There are four G130M observations of 
this sightline, with central wavelengths of ISISA, 1309A, 
1300A, and 129lA respectively, yielding a total integration 
time of 4401 s. We generate a series of bins approximately 21 
pixels in size (i.e. three resels, or 0.21 A), and within each bin 
we average together the fluxes from each pixel of each spec¬ 
trum that falls within the bin. The merged G130M nuclear 
spectrum is shown in Figure 7. 

This spectrum contains many well-known lines typi¬ 
cal for photoionized or collisionally excited plasmas with 
Tf, ~ 10^ — 10® K , many of which have optical counterparts 


that have previously been studied with HST (e.g. [Ford et 
al.||1994| [Harms et al.|[l994[ [Macchetto et al.[|1997| [Walsh 


et al. 


20131. Of particular note, the emission lines (partic¬ 


ularly NV A1240 doublet and the Gii A1336 triplet) are 
extremely broad and C ll appears to have a blueshifted ab¬ 
sorption component as well. Our interest in this paper is 
primarily with the [Fexxi] line, which uniquely traces the 
10^ K gas (note that [Fexix] falls outside the spectral cov¬ 
erage of these G130M observations). 

Using the merged spectrum from Figure 7, we repeat 
the Gaussian fitting procedure from Section 2 at the loca¬ 
tion of [Fexxi]. We find vq = —10 km/s, a = 290 km/s, and 



Figure 6. Estimated emissivities (erg cm® s“®) for the sum of 
the emission from the Civ A1549 doublet (blue) and the sum of 
the Hell A1640 multiplet (red), assuming Solar abundance and 
CIE. The predicted ratio of these lines at 10® K is 110; to get the 
observed ratio of ~ 1.0 a significant contribution from hotter gas 
is required. The ratio from cooling flow models is still too large 
by nearly an order of magnitude, but thermal conduction models 
seem to be able to reproduce the observed ratio ([Sparks et al.[ 
\ 2012 \ . 


F = 1.4 X 10“^® erg s“^ cm“®. The best-fit profile is shown 
in Figure 7 (lower right panel). This is not a significant de¬ 
tection of [Fexxi] (S/N = 0.7). The 90% upper limit on the 
integrated flux is 4.15 x lO”'^® erg s“^ cm“^, corresponding 
to a line luminosity of L < 1.4 x 10®^ erg s“^ within this 
sightline. 

In Figure 8, we repeat the emission measure analysis 
which we performed above for the filament, using the same 
pressure of P = 1.91Pe = 4x 10® cm“® which we used for the 
filament region. The l.”25 radius of the COS aperture cor¬ 
responds to 98 pc here, which is roughly the same size as the 
Bondi radius for M87 ( |Di Matteo et al.j2003| ), although the 
sightline is offset from the black hole so the plasma within 
the Bondi radius only covers a portion of our aperture. The 
pressure within the Bondi radius increases above our fiducial 
value, but the emission measure only rises logarithmically 
towards the center ( Russell et al.|2015 |, so the contribution 
from the Bondi flow to [Fe xxi] in this region is not likely to 
be dominant. 

We find that the 90% upper limit on the minimum emis¬ 
sion measure is 7.4 x 10®® cm“®. The electron column is 
Ne < 4.7 X 10®® cm“®, the effective path length is I < 8.2 
kpc, and the hot gas mass is M < 1.5 x 10^Mq within our 
pencil-beam aperture. 

There is also an X-ray measurement of the emission 
measure of 1 keV gas in this region, from the 5-temperature 
spectral modeling for this region performed by [Werner "et] 


al. (20101. From their high-resolution Chandra map, they 


estimate the emission measure of 1 keV gas from the nucleus 
of M87 to be around 6 x 10®® cm“® arcsec”®, assuming a 
distance of 16.1 Mpc and a metallicity of 1.3^0. Scaling this 
measurement to an assumed Solar metallicity, we predict 
an [Fexxi] A1354.1 flux of 7 x (Z/Zq) erg s“® cm“® 
based on these X-ray observations. This is consistent with 
our observations. 
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Figure 7. (top) Average spectrum of four COS G130M observations of the nucleus of M87. The spectrum is shown in black and the 
1(7 error per bin is shown in green. The spectrum has been binned to a resolution of 21 pixels (approximately three resels). Prominent 
airglow lines are indicated in light blue. The Lyman-c« aigrlow feature is blended with the extremely broad Lyman-a emission line from 
the nucleus of M87 and we truncate these lines for clarity. Other lines are also identified and marked; some of them are also extremely 
broad. The lower left panel shows the spectrum within 2500 km/s of [Fexxi], as well as our linear fit to the continuum. There is no 
evidence for geocoronal OI contamination in this region. The lower right panel shows the background-subtracted spectrum within 1000 
km/s of [Fexxi] as well as the best-fit Gaussian profile, which corresponds to a 0.7(7 detection. The 90% upper limit on the integrated 
[Fexxi] line flux is 4.15 X 10“^® erg s“^ cm“^. 


4 NGC 4696 NUCLEUS 


We also examine a sightline towards the center of NGC 4696. 
This is the brightest central galaxy of the Centaurus cluster, 
which is a classic cool-core system: the temperature of the 
intracluster medium decreases from a baseline of 3-4 keV 


down to a bit below 1 keV in the center (Fukazawa et al. 


1994| [Allen & Fabian||19941 [Allen et al. 

20011 IMolendi et 

al.|2002| Fabian et al.|2005 [Graham et a 

i.|2006| Sanders & 


Fabian|2006 Sanders et al.|2008 Panagoulia et al.|2013 1. As 
with M87, the X-ray spectra in the cool core prefer multi¬ 
temperature fits, and NGC 4696 also hosts extensive fila¬ 
mentary structure which radiates in optical emission lines 


(Fabian et al. 1982 Crawford et al. 20051. More recently, 
[Fex] A6374 emission has also been detected from t his galaxy 
( Canning et al.|[2011 |, which traces gas at 10® K. [Mittal et 
( |2011[ ) have also detected [C ll] A158/rm, [O l] A63/rm, 
and [N ll] A122/rm emission which is largely coincident with 
the filamentary structure, suggesting the presence of 10^ K 
gas as well. The sightline towards this galaxy also covers its 
nucleus, which is an FR I radio source. 

The four datasets covering this sightline were taken in 


a similar configuration as the observation of the nucleus of 
M87, and have a total integration time of 5217 s. We average 
the four spectra and bin the resulting spectrum to 56 pixels 
(eight resels) per bin, in the same fashion as the nucleus of 
M87. NGC 4696 has a redshift of g = 0.009867 (u = 2958±15 
km/s; de Vaucouleurs et al.|[l9^ l, which we use to correct 
the spectra to the rest frame. Using the updated Planck 
cosmology ( [Planck Collaboration et al.||20i^ this gives a 
luminosity distance of 44.0 Mpc and implies a radius for the 
COS aperture corresponding to about 260 pc for this galaxy. 


The resulting merged spectrum is shown in Figure 9. 
This spectrum has a much lower continuum than the spec¬ 
trum of the nucleus of M87; this is partially due to the higher 
distance of NGC 4696, and partially intrinsic. This observa¬ 
tion is also much more heavily contaminated by airglow than 
the M87 observations: in addition to prominent Hi emission 
and strong emission from both Ol features, the Nl triplet 
at 1200A is also visible. Fortunately, the higher redshift of 
this galaxy means that none of these airglow lines fall near 
[Fe xxi]. 

We repeat the analysis we performed for M87. We 
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Figure 8. The 90% upper limit on the pencil-beam emission 
measure f neUpdl from the nuclear sightline towards M87 inferred 
from our limits on the measured [Fexxi] flux. We also compute 
limits on the corresponding average electron density (upper X- 
axis), electron column density (dotted lines), and path length 
(dashed lines) for the hot gas, assuming a constant t otal pressure 
of P = 1 .91Pe = 4 X fO® cm“® K in the sightline ( [Churazov et| 
|al.||2008| . This curve is sharply peaked at 10^'®® K, and if the 
[Fe xxi]-emitting gas has this temperature the limits approach 
physically plausible parameters for the hot gas in the sightline. 
The star indicates the approximate emission measure of 1 keV gas 
estimated from the X-ray observations of [Werner et al.| l |2010| | , 
which we have scaled to an assumed Solar metallicity. This is 
consistent with our limits. The physical values in this plot are all 
computed assuming Solar abundance, and scale with metallicity 
as discussed in the text. 


fit a linear continuum to the region around the line. Our 
best-fit Gaussian has uo = 0 km/s, cr = 300 km/s, and 
F = 2.0 X 10“^^ erg s“^ cm~^, but it is only significant at 
0.2(7. The more relevant quantity is the 90% upper limit on 
the integrated flux, which is 2.2 x 10“'^® erg s“^ cm“^ (cor¬ 
responding to a line luminosity L < 5.1 x 10®^ erg s“^). We 
are somewhat sensitive to an accretion flow in this galaxy 
as well, but NGC 4696 has a less massive black hole than 
M87 and its Bondi radius fills a much smaller fraction of our 
aperture than in the M87 observation. 

In Figure 10, we repeat the emission measure analysis 
which we performed above for M87. Here we assume a pres¬ 
sure of 2 X 10® cm“® K based on the deprojected density and 
temperature profiles at 1 kpc measured by [Graham et al.j 
(20061. We could not find estimates in the literature against 
which to compare the emission measure within our aperture, 
but we find similar limits on the [Fe xxi]-emitting column in 
NGC 4696 as compared to M87. The upper limit on the 
minimum emission measure is 3.1 x 10^® cm“®. The electron 
column is Ne < 4.0 x 10^^ cm“^, the effective path length 
is / < 14.1 kpc, and the hot gas mass is M < 8.3 x 10 ®Mq 
within our pencil-beam aperture. 


5 POTENTIAL SOURCES OF 
CONTAMINATION 

In this section, we briefly explore other possible explana¬ 
tions for the positive residuals seen in our spectra towards 


M87 and NGC 4696. As discussed in Section 2 and in the 
Appendix, at the redshift of M87, [Fe xxi] falls near the geo¬ 
coronal OI A1356 lines, although in most of these observa¬ 
tions the lines do not appear to contribute significantly. 

It is also important to understand the continuum in 
these spectra, since the continuum emits much more flux 
than the forbidden lines we consider. M87 is one of the classic 


elliptical galaxies with an upturn in the FUV (Bertola et 


al. 19801. The exact cause of this ”FUV upturn” is still 


unclear, but the likeliest possibility seems to be some type 
of hot evolved star, either isolated or in a binary system (see 
the review by 0’Connell||1999 more discussion). Importantly, 
the FUV upturn, when it exists, appears not to have any 
prominent emission lines (Ferguson et al.|1991[ Brown et al. 


1995 Brown et al. [[19117 ( and might even have a few broad 


absorption features ( |Bica et al.|1996 (. We can therefore treat 
the continuum as being featureless, to first order, in the 
region around [Fexxi] and [Fexix[. NGC 4696 has a weaker 
continuum, but it likely has a similar origin. 

Stellar coronae of some late-type stars also show emis¬ 
sion in [Fexxi] A1354.1, and in principle these may con¬ 
tribute to the flux in the lines of interest. In the Capella bi¬ 
nary, each star ha s an [Fexxi] A1354.1 luminosity of around 
4 X 10^® erg s“^ (Linsky et al. 19981, so our upper limits 
correspond to ~ 10^'^ — 10^^ Capella-type giant stars in each 
sightline, which is orders of magnitude above the expecta¬ 
tion. We therefore do not expect stellar coronae in active 
binaries to contribute a major fraction of the [Fexxi] from 
these sightlines. 

Observations should be able to separate the contribu¬ 
tions from stars and gas. The stars in the nuclei of M87 have 
a velocity dispersion of 400-450 km/s ( Gebhardt et al.|2011 ( 
and in NGC 4696 they have a velocity dispersion of 250 km/s 
(HyperLEDA; [Makarov et al.|[2014 1 respectively. Our best- 
fit Gaussian profiles are broad, but in the M87 nucleus we 
do not see this level of broadening, and deeper observations 
should be able to distinguish a stellar velocity dispersion 
from the smaller velocities expected from Iron ions. 

Finally, in nearby stellar coronal spectra that show the 
[Fexxi] line, the nearby Cl A1354.3 line can often be an 
issue. Based on non-detections of the other C l lines in these 
spectra it is unlikely that Cl A1354.3 is significant here, 
and if it were present we would again expect it to have the 
velocity dispersion of the stars. 

On the other hand, for an early-type galaxy without an 
FUV upturn or any 10^ K gas, extragalactic observations 
which target the [Fe xxi] emission from coronally-active late- 
type stars are conceivable. Such observations would be a 
unique constraint on the properties of these stars in elliptical 
galaxies. 


6 OTHER SIGHTLINES 


Here we briefly consider additional sightlines towards cool¬ 
ing flows or filaments seen in massive galaxy clusters. These 
clusters are all hotter than 1 keV, so they are not expected 
to host significant reservoirs of 10^ K gas, but our sight¬ 
lines do probe regions of multiphase cooling so some 10^ K 
gas likely exists. We found publicly available COS observa¬ 
tions towards A1795 {z = 0.0619; McDonald et al. 20141, 
Perseus {z = 0.0176, Churazov et al.||2003 (, Phoenix (SPT- 
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Figure 9. (top) Average spectrum of four COS G130M observations of the nucleus of NGC 4696. The spectrum is shown in black 
and the la error per bin is shown in green. The spectrum has been binned to a resolution of 56 pixels (approximately seven resels). 
Prominent airglow lines are indicated in light blue, and truncated for clarity; airglow contamination is much more significant in this 
spectrum than in the observations of M87. The lower left panel shows the spectrum within 2500 km/s of [Fexxi], as well as our linear fit 
to the continuum. There is weak geocoronal OI emission in this region. The lower right panel shows the background-subtracted spectrum 
within 1000 km/s of [Fexxi] as well as the best-fit Gaussian profile, which corresponds to a O.lcr detection. The 90% upper limit on the 
integrated [Fexxi] line flux is 2.2 X 10“^® erg s“^ cm“^. 


CL J2344-4243; 2 = 0.596; McDonald et al.]|2012 |, and Zw 
3146 {z = 0.2906, Allen et al.||1992 i. We also consider the 
COS observations of a shock front in HCG 92 {z = 0.021; 


White et al. 19991, which is likely cooler than 10^ K, but 
may contain a hotter component at 10^ K. 

The sightline towards A1795 targets a filament pro¬ 
jected about 30 kpc from the central galaxy. According to 
deep X-ray observations (Ehlert et al. 20151, the filament 
can be fit with a multi-temperature model, but the hotter 
component is around 4 keV and the cooler component is 
about 2 keV, which is not quite ideal for producing these 
FUV lines. However, the COS G140L spectrum shows a 
weak positive residual at [Fexxi], significant at O.lcr, with 
a 90% upper limit on the integrated flux of 1.3 x 10“^® erg 
s“^ cm“^. The deprojected electron pressure at the location 


of this sightline is around 8 x 10® cm ® K (Tamura et al. 


2001 Fttori et al. 20021; with this assumption and an as¬ 


sumed Solar metallicity our upper limit on the effective path 
length of gas at 10^'®® K is / < 23 kpc and the upper limit 
on the electron column density is W < 5 x 10^^ cm“^. 

In the A1795 spectrum there is also a weak detection 
of Ovi A1038 claimed by McDonald et al. (20141. They 


place a upper limit of 90 km/s on the width of the line, 
which would imply that the line is not significantly spatially 
broadened (recall from section 1.1 that aperture-filling emis¬ 
sion observed with G140L can be expected to have cr ~ 650 
km/s). At the redshift of A1795, the COS aperture subtends 
a diameter of about 3 kpc, so if the filament is significantly 
narrower than this, spatial broadening might not be impor¬ 
tant. On the other hand, the thermal velocity of O ions at 
10® K is only about 10 km/s, so there remains the possibil¬ 
ity of up to 90 km/s of turbulent broadening and/or spatial 
broadening in this feature. 


In Perseus we examine a deep observation of the so- 
called outer “horseshoe” filament projected about 20 kpc 
from the center of NGC 1275. The temperature of the intr¬ 
acluster medium at this location is around 3 keV, much too 
hot to expect emission from [Fexxi], and indeed we see no 
trace of this line in the combined G140L spectrum. The 90% 
upper limit on the flux is 6.2 x 10“^^ erg s“^ cm~^. Com¬ 
bining th is with the estimated e lectron pressure of 1.4 x 10® 
cm“® K (Churazov et al. 20031 at this location, we derive 
an upper limit on the effective path length of gas at 10^'®® 
K of / < 3 kpc and an upper limit on the electron column 
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Figure 10. The 90% upper limit on the pencil-beam emission 
measure f neUpdl from the nuclear sightline towards NGC 4696 
inferred from our limits on the measured [Fexxi] flux. We also 
compute limits on the corresponding average electron density 
(upper X-axis), electron column density (dotted lines), and path 
length (dashed lines) for the hot gas, assuming a constant total 
pressure of P = 1.91Pe = 2 X 10® cm“® K in the sightline ( |Gra-| 
|ham et al.|2006| . This curve is sharply peaked at 10^ ®® K, and if 
the [Fe xxi]-emitting gas has this temperature the limits approach 
physically plausible parameters for the hot gas in the sightline. 
The physical values in this plot are all computed assuming Solar 
abundance, and scale with metallicity as discussed in the text. 


hand, implies a [Fexxi] A1354.1 limit of around 6 x 10“^® 
erg cm“^. Our COS limits are much deeper than this, 
so if this redshifted line fell in the G160M bandpass COS 
would be able to improve considerably the constraint on 10^ 
K gas in the cool core of Phoenix. 

Zw 3146 also has a moderate redshift and its G140L 
spectrnm has good sensitivity for [Fexxi], [Fexix], and 
[Fexxillj. This clnster has a core temperature of around 3 
keV ( [Diehl et al.||2008 l and a cooling flow nearly as strong 
as Phoenix ( Edge et al.[|1994 1, and the COS sightline passes 
through this cooling flow. These three lines are not de¬ 
tected in this spectrum; the respective 90% upper limits are 
2.4 X 10“^® erg s“^ cm“^, 1.2 x 10“^® erg s“^ cm“®, and 
0.8 X 10“^® erg s“^ cm“®. Again we do not have a good 
constraint on the hot gas pressure within this sightline, so 
we simply report the upper limits on the flux here. 

Finally, HCG 92 (Stefan’s Quintet) has a diffuse intra¬ 
group medium with kT « 0.5 keV ( Trinchieri et al.]|2003 l, 
and a large X-ray bright shock near NGC 7318b which is 
the target of the COS observations. Temperature measure¬ 
ments are often difficult in shocked plasma, but [Trinchieri etj 
E] ( [2005P infer a multi-temperature medium with the hot¬ 
ter component around 1 keV while O’Sullivan et al. (20091 


density of A^e < 1 x 10 cn 
limits. 

The Phoenix cluster is located at much higher redshift 
than the other targets we have considered so far. This allows 
us to constrain some of the other lines in Figure 1, which lie 
at wavelengths too short for COS to observe with good sen¬ 
sitivity in the local Universe. As long as our source fills the 
COS aperture, the 1/R^ decrease in the flux is matched by 
a 1/7?® increase in area, so our sensitivity remains roughly 
constant (except for a small decrease stemming from the 
different redshift dependences of the angular diameter dis¬ 
tance and the luminosity distance), allowing us to derive 
useful constraints even at moderate redshifts. 

The temperature of this supermassive cluster is around 


hnd a single-temperature model at about 0.6 keV is suffi¬ 
cient. We combine all the available COS data, which targets 
multiple sightlines along the X-ray shock front, and hnd no 
evidence for any of our high-excitation FUV lines. Unfor¬ 
tunately [Fexxi] lies in the gap between the G130M and 
the G160M observations, but we obtain 90% upper limits 
for [Fexix] and [Fexxiii] of 3.7 x 10“^® erg s“^ cm“® and 
which are both fairly tight 7.4 x 10“^® erg s“^ cm“® respectively. Given the uncertain 


13 keV (McDonald et al. 20121, but the COS observations 
target the powerful (M ~ 38OOM0 yr“^; McDonald et al. 


properties of the gas in the shock, the pressure is not well- 
constrained, and we simply report the upper limits on the 
hux. 

In Table 3 we collect many of the constraints reported 
in this work. For each FUV line we examine, we list the 
constraint on the flux and the associated pencil-beam emis¬ 
sion measure.The latter quantity depends on the metallicity 
and temperature of the gas; we adopt the same convention 
as used above in our emission measure analyses and assume 
Solar abundance and that each species is at the temperature 
corresponding to its peak emissivity. Defined in this way, the 
quoted emission measure can be derived from the flux with¬ 
out any information other than an assumed distance to the 
target, which we derive from the quoted redshift in every 
case except for M87, where we use a distance of 16.7 Mpc 


20121 cooling flow. X-ray observations of the cooling flow 


show a multiphase medium, but they do not detect any gas 
below 1.8 keV, with upper limits on the cooling rate below 
this temperature on the order of hundreds of Mq yr“^ Tozzi 


[et al.[ ( [2015 1. The G160M spectrum covers [Fexviii] and 
[Fe xxill], but neither line is detected. Our 90% upper limits 
on the flux are 1.7 x 10~^® erg s“^ cm“® and 1.9 x 10“^® 
erg s“^ cm~® respectively. These limits are comparable to 
the X-ray limits from the Tozzi et al. ( 2015[ ) RGS spectrum. 
While a detailed comparison would be outside of the scope 
of this paper, the ratio of Fexxi A12.16 to [Fexxi] A1079 
is about 50 at 10^ K, so our upper limit corresponds to an 
upper limit of 1 x 10“^"^ erg s“^ cm~® for the X-ray line. 


From Fig. 6 of Tozzi et al. (20151, the integrated flux of this 
line is sightly under this limit (we estimate 4 x 10“^® erg 
®). The non-detection of Fexxi A9.48, on the other 


(Blakeslee et al. 20091. In the systems where we have an 
estimate of the hot gas pressure, we also list the derived 
constraint on the effective path length of hot gas and the 
free electron column density. 

Assuming we are at the the peak of the emissivity for 
each line, the emission measure / neUpdl is the fundamental 
parameter and is linearly proportional to flux, so deeper 
observations will improve these limits proportionally to the 
square root of the exposure time. For many of the objects 
discussed in this section, such as Phoenix and HGG 92, the 
integration time would have to increase by factors of several 
in order to obtain emission measure limits as good as some 
of the other targets in our sample; a better approach would 
be to focus on [Fe xxi] for these targets since it has a higher 
emissivity than the other lines. 

On the other hand, the derived column density and ef- 
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Table 2. 90% Upper Limits Reported in this Work 


Target 

Line 

Integrated Flux 

(IQ-ie gj,g g-i cm“^) 

j TleTlpdl 
(10^^ cm-^) 

1 

(kpc) 

Ne 

(10^1 cm-2) 

M87 filament (night-only) 

[Fe xxi] 

1.7 

3.2 

3.5 

2 

M87 filament (day+night)* 

[Fe xxi] 

0.8 ± 0.4 

1.4 ±0.6 

1.5 ±0.7 

9±4 

M87 filament (day+night) 

[Fe xix] 

1.5 

5.9 

3.8 

3 

M87 nucleus 

[Fe xxi] 

4.2 

8.0 

8.2 

5 

NGC 4696 nucleus 

[Fe xxi] 

2.2 

4.1 

14.1 

40 

A1795 filament 

[Fe xxi] 

1.3 

2.9 

23 

50 

Perseus filament 

[Fe xxi] 

0.6 

1.2 

3 

10 

Phoenix nucleus 

[Fe xviii] 

1.7 

23 

- 

- 

Phoenix nucleus 

[Fe xxiii] 

1.9 

310 

- 

- 

Zw 3146 nucleus 

[Fe xxi] 

2.4 

12 

- 

- 

Zw 3146 nucleus 

[Fe xix] 

1.2 

12 

- 

- 

Zw 3146 nucleus 

[Fe xxiii] 

0.8 

55 

- 

- 

HCG 92 shock 

[Fe xix] 

3.7 

13 

- 

- 

HCG 92 shock 

[Fe xxiii] 

7.4 

27 

- 

- 


90% upper limits reported in this work. For each line we constrain along each sightline, we present the 90% upper limit 
on the integrated line flux as well as the upper limit on the pencil-beam emission measure. For the latter quantity, 
we assume Solar abundance and that the gas is at the temperature corresponding to the peak of the emissivity for 
that line. In systems where we have independent constraints on the gas pressure (from e.g. X-ray observations) we 
combine this information with the emission measure constraints to infer upper limits on the effective path length I 
and the free electron column density Ne of the hot gas as well, using equations 1 and 2. *For the M87 filament, we 
report both a 90% upper limit based on a night-only spectrum as well a 2.2 (t detection of the line which includes 
observations taken during orbital daytime. For the latter values, Itr uncertainties are also listed. 


fective path length also depend on the pressure of the hot 
gas (see eqs. 1 and 2). The latter quantity in particular de¬ 
pends on the inverse square of the pressure, so even weak 
constraints on the emission measure can translate to tight 
constraints on the effective path length if the ambient pres¬ 
sure is high (i.e. in the cores of massive galaxy clusters). 

Of course, the best target in this list is the filament 
region in M87, where we have a tentative 2a detection of 
[Fexxi]. In the next section, we show that modest follow¬ 
up observations should readily improve these statistics and 
allow us to measure line kinematics. We also consider the 
nuclei of M87 and NGC 4696 as well. 


7 UPDATED HST PIPELINE 

After the bulk of this analysis had been completed, the new 
HST Spectroscopic Legacy Archiv^was released, which of¬ 
fers science-grade co-added spectra for galaxies and galaxy 
clusters observed with COS. In general these combined spec¬ 
tra agree well with the co-added spectra produced in this 
paper, which are based on the calcos pipeline. However, 
the HSLA uses a different (and less conservative) method 
of estimating the errors, resulting in significantly smaller 
flux uncertainties than the errors produced by the calcos 
pipeline. We therefore repeated the searches for forbidden 
Iron lines for each spectrum using these new HSLA spectra, 
and we report the results in this section. 

For the filament region in M87, the HSLA full spectrum 
yields a positive residual at the expected location of [Fe xxi] 
with Vo = —560 ± 40 km/s, a = SSO^go km/s, and an inte¬ 
grated flux of 9.6±1.8x 10“'^^ erg s“^ cm~^. This is generally 

^ http://archive.stsci.edu/hst/sp6ctral_legacy/ 


similar to the best-fit values in section 2, although the veloc¬ 
ity is blueshifted by 300 km/s compared to the previous fit; 
this likely reflects uncertainty in continuum fitting due to 
the low S/N of the spectrum. The formal significance of the 
detection is now 5.3a (compared to 2.2a from our estimate 
in section 2), although this should be reduced somewhat due 
to uncertainties in the continuum. 

We also estimated the line properties in the night-only 
spectrum using the HSLA error estimate. The night-only 
spectrum is not provided by the HSLA but we were able to 
reproduce their algorithm for computing the errors, which 
we then applied to our co-added night-only spectrum. We 
find Vo = —4401^20 km/s, a — 390/(f jq km/s, and an inte¬ 
grated flux of h.blj'g 10“^^ erg s“^ cm“^. These values are 
again reasonably consistent with the values from our previ¬ 
ous analysis, albeit with higher total significance for the line 
(S/N = 2.2 for this analysis, as compared to 1.0 in section 
2 ). 

Under the new error estimate, the [Fexix] line in the 
M87 filament has vo = — 130429o km/s, a = 3504^20 km/s, 
and an integrated flux of F = x 10“^^ erg s“^ cm~^. 

Again, these values are consistent with our previous esti¬ 
mates, but the S/N has increased from 0.3 to 0.8. 

For the nucleus of M87, the feature at the expected 
location of [Fexxi] now has vq = O/iiao km/s, a = 270^^^ 
km/s, and an integrated flux of 1.04 ± 1.44 x 10“^® erg s“^ 
cm“^. These parameters are very similar to their previous 
best-fit values; the uncertainties are slightly decreased but 
the significance of the detection remains constant at 0.7 ct. 

For NGC 4696, we find for [Fexxi] uo = 0 ± 130 km/s, 
a = 28OI230 km/s, and F = 0 . 8 /(g | x 10 “^^ erg s“^ cm~^. 
This is consistent with our previous result for the line, but 
the S/N has actually decreased from 0.3 to 0.1. The 90% 
upper limit on the flux is 1.28 x 10“^® erg s“^ cm“^, which 
is slightly weaker than the previous constraint. 
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The upper limits in section 6 remain upper limits, but in 
general these limits are tighter using the HSLA spectra. The 
significance of the positive residual seen towards the filament 
in A1795 increases from 0.7a to O.Scr, and the 90% upper 
limit on the integrated flux decreases to 7.4 x 10“^^ x 10“^^ 
erg s“^ cm~^. In the “horseshoe” filament of Perseus, our 
90% upper limit also tightens, to 2.6 x 10“^^ erg s“^ cm“^. 
In the Phoenix cluster, our limits on [Fe xvill] and [Fe xxill] 
become 1.2 x 10“^® erg s“^ cm“^ and 1.3 x 10“^® erg s“^ 
cm“^ respectively, which are moderate improvements over 
the values in section 6. For ZwCl 3146, the 90% upper lim¬ 
its on [Fexxi], [Fexix], and [Fexxiii] respectively are now 
1.1 X 10“^® erg s“^ cm“^, 0.9 x 10“^® erg s“^ cm“^, and 
0.8 X 10“^® erg s“'^ cm“^. The first two values are moderate 
improvements over the limits in section 6, while the third 
is unchanged. Finally, in HCG 92, the HSLP data do not 
include the G130M data, so we are unable to update the 
limits for this sightline using the new error method. 

In general, using the new HSLP method of calculating 
errors yields similar results to the calcos pipeline errors, but 
with higher inferred significances for detections and tighter 
constraints on non-detections. To remain conservative we re¬ 
port results based on the latter errors in this work, but it 
is straightforward to replace the values in Table 2 and else¬ 
where with the results based on the HSLP errors if desired. 


8 DETECTABILITY IN FUTURE 

OBSERVATIONS 

The required integration time depends on the width of the 
line, and we consider four example cases for the velocity 
dispersion: cr = 50 km/s (the thermal velocity of Iron ions), 
a = 150 km/s (30% of the sound speed for 10^ K gas), 
a = 500 km/s (the sound speed for 10^ K gas), and a is 
equal to the velocity dispersion of the stars in the galaxy 
(roughly 450 km/s for M87 and roughly 250 km/s for NGG 
4696). This yields a FWHM for [Fexxi] of 0.53A, 1.59A, 
5.31A, and 4.78/2. 66A (for M87/NGG4696) respectively for 
the four fiducial values of a. 

Using our estimated flux of 7.8 x 10~^^ erg s“^ cm“^ 
for [Fexxi] from the filament in M87, the times required 
to achieve an integrated S/N of 5 would be about 6 ks, 
12 ks, 33 ks, and 30 ks respectively, or approximately 3, 
6, 11, and 10 orbits of integration. This is comparable to 
the existing integration time with the G140L grating, but 
the medium-resolution gratings are more sensitive and the 
spatial broadening will be less severe with these gratings as 
well. 

In the M87 nucleus, our predicted flux is similar (7 x 
10“^^ erg s“^ cm“^ for Solar metallicity; see Section 3). 
However, the high continuum in this region introduces noise 
as well. The continuum is about 1.3 x 10“^® erg s“^ cm“^ 
A“^, and it dominates the noise budget. In order to detect 
the line with an integrated S/N of 5, for the four fiducial line 
widths we need S/N values of about 17, 29, 54, and 48 per 
resel for the continuum. This requires about 40 ks, 110 ks, 
390 ks, and 310 ks respectively. The latter three integration 
times can only be obtained through a Large GO program. 
It would probably be much easier to observe a sightline a 
few arcseconds away from the center of M87, where the con¬ 
tinuum is much lower. If the line flux remained constant 


but there were no continuum, the required integration time 
would be reduced by a factor of 10. 

In the NGG 4696 nucleus, the continuum is less of an 
issue, but the line flux is very low. If we use the best-fit 
integrated flux of 2 x 10“^^ erg s“^ cm“^ for [Fexxi] in this 
sightline, then the required integration times to achieve a 
S/N = 5 are 50 ks, 130 ks, 410 ks, and 220 ks respectively. 

The velocity resolution per resel is around 15 km/s for 
G130M, and this is also about the la absolute calibration 
uncertainty for the instrument. However, depending on the 
width of the line, some binning may be necessary to get a 
S/N of 1 per bin if the integrated S/N is 5. Accounting for 
this binning, for our four fiducial velocity dispersions, we 
expect velocity resolution of about 20, 30, 60, and 40/60 
km/s respectively. 

If the emission is extended, then this velocity resolution 
will be degraded further by spatial broadening. In section 
1.1 we estimated this effect to be about 81 km/s for uni¬ 
formly extended emission. Extended emission can be iden¬ 
tified from the 2D cross-dispersion axis versus dispersion 
axis image, and deconvolution can be attempted if neces¬ 
sary, but a limiting resolution of 81 km/s is still just 15% of 
the sound speed in 10^ K gas. A measurement of turbulent 
motion with ~80 km/s resolution would be better than any 
previous measurements. 

It is worth emphasizing that FUV spectrographs have 
much better sensitivity than their X-ray counterparts, which 
enables the study of these faint lines. For example, the Re¬ 
flection Grating Spectrometer on XMM-Newton has an ef¬ 
fective area of up to ~ 150 cm^ and the Soft X-ray Spectrom¬ 
eter on Hitomi had ~ 200 cm^, while GOS has more than 
2000 cm^ with the medium-resolution gratings. Proposed 
future UV spacecraft such as the World Space Observatory 
) or the Public Telescop^ will retain 
or improve upon the spectral resolution of GOS, but these 
missions do lack the large collecting area of HST. 

Finally, as we mentioned in the introduction, these lines 
can be observed in absorption as well. The absorption cross- 
sections are low, typically of order or, but with absorption¬ 
line observations much more distant sources can be analyzed 
as long as they have a bright background quasar projected 
behind them. For [Fexxi], the optical depth is 1% for an 
electron column of 8 x 10^^ cm“^ (assuming Solar metallicity 
and a temperature of 10^'®® K). This can be achieved with 
an electron density of 0.2 cm“® over a path length of 10 kpc, 
which is obtainable in the center of some galaxy clusters, for 
example. 


9 CONCLUSIONS 

In this paper we searched for [Fexxi] A1354.1 and [Fexix] 
A1118.1 emission in publicly available COS FUV spectra to¬ 
wards M87 and NGG 4696. These spectra are very complex 
and we do not attempt to model them in full, instead focus¬ 
ing our analysis on the region around [Fexxi] and [Fexix]. 
Towards a filament in M87, we report a tentative (2.2cr) de¬ 
tection of [Fe xxi] with an integrated flux of 7.8±3.6 x 10“^^ 
erg s“^ cm“^. This line is located very close to a nearby O l 

^ http://www.publictelescope.org/ 


Shustov et al. 2009 
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clei of M87 and NGC 4696, with values of 4.5 x 10 


airglow line, and while this airglow line seems to be neg¬ 
ligible in these observations, we also report a 90% upper 
limit based only on the data taken during orbital night. We 
also report 90% upper limits on [Fexxi] towards the nu- 

erg 

erg s^" cm^* respectively, corre¬ 
sponding to limits on the line luminosities of 1.4 x 10®^ erg 
s“^, and 5.1 x 10®^ erg s”'^. 

We perform emission measure analysis on each pencil- 
beam column as well, inferring constraints on the electron 
column and the characteristic path length of the hot gas, 
assuming a temperature of 10^'°® K (the temperature at 
which [Fexxi] emissivity peaks). For the filament sightline, 
we infer a hot gas electron column of 9 ± 4 x 10^° cm“^ and 
an effective column length of 1.5 ± 0.7 kpc. For the nuclei 
of M87 and NGC 4696, our respective 90% upper limits on 
the column density are 7.4 x 10^° cm~^ and 4.0 x 10^^ cm“^ 
and our respective 90% upper limits on the effective path 
length are 8.2 kpc and 14.1 kpc. 

In both the filament sightline and the nuclear sightline 
we are able to compare our constraints to X-ray measure¬ 
ments of the emission measure of 1 keV gas by [Werner et al.| 
(20101 and Werner et al. (20131. In both cases our UV con¬ 
straints are fully consistent with the X-ray measurements. 

We also examine several other sightlines towards galaxy 
clusters. These additional sightlines probe multiphase me¬ 
dia, but there are not clear predictions of massive reservoirs 
of 10^ K gas along these sightlines. We do not detect any 
of these FUV lines in these spectra, although some of our 
upper limits are competitive with the M87 and NGC 4696 
spectra, or with existing X-ray measurements. Our limits 
can be used to constrain the differential emission measure 
of the multiphase plasma in some cases. 

Follow-up COS observations of moderate duration 
should be able to improve the situation dramatically, as we 
showed in Section 8. In particular, it is especially important 
to obtain medium-resolution spectra of the M87 filament, in 
order to confirm the weak [Fexxi] feature we identified and 
to measure the shape of the line. 

With future instrumentation, observations targeting 
[Fexxi] will be able to expand to additional sightlines. Mul¬ 
tiple sightlines for a single object could be extremely valu¬ 
able for mapping out the velocity structure of the intraclus¬ 
ter and intragroup media, as well as quantifying the density 
fluctuations ( Churazov et al.||2012 Zhuravleva et al^|2014 |. 
Additionally, constraining the other FUV lines in Figure 1 
along with [Fexxi] would yield precise constraints on the 
distribution of temperatures in the hot gaseous halo. 
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HST, and so the emission from this doublet looking towards 
zenith is generally negligible ( Syphers et al.|[2012 |). As the 
altitude increases, the zenith-looking ratio of Ol A1356 to 
OI A1304 gets lower and lower, since the latter have oscilla¬ 
tor strengths which are orders of magnitude higher than the 
former. Some O i A1304 emission produced below the orbit of 
HST can therefore scatter into the line of sight, while most 
nightglow OI A1356 emission does not reach the detector 
( Gerard et al.^loT? Chakrabarti et al. 19841. This explains 
the relative nightglow line strengths listed in the Gycle 23 
COS instrument handbook (table 7.4), in which Ol A1304 
is an order of magnitude brighter than O l A1356 at orbital 
night. 

In Figure 2, the O l A1304 triplet has nearly disap¬ 
peared, due to the removal of the dayglow emission. The 
emission from the Ol A1356 should be even weaker than 
the tiny residual at 1304A, probably by an order of magni¬ 
tude or more. Consequently, we can feel confident that the 
spectrum in Figure 2 is essentially clean of any measurable 
O lA1356 geocoronal emission. 

The analysis of dataset LC7T02010 was more compli¬ 
cated, so we do not include results from this observation 
in the merged spectrum in Figure 2. Three of the ten ob¬ 
servations within this dataset (Ic7t02h0q, Ic7t02hiq, and 
Ic7t02hmq) show anomalously high O l A1354 airglow emis¬ 
sion, which remains strong during orbital night even while 
011304 decreases. This behavior has been observed by other 
missions as well (e.g. DeMajistre et al.|2005 1 during periods 
of high geomagnetic activity. We tested the effect of combin¬ 
ing the other seven observations in dataset LC7T02010 with 
the observations in Figure 2, but the restriction to orbital 
night removes most of the data, so that only 883 seconds of 
integration time remain to be added to the existing 7858 s 
of night-only data. The results are therefore nearly identical 
to those shown in Figure 2. 


APPENDIX A: FILTERING O I 1356 AIRGLOW 
EMISSION FROM THE G140L SPECTRA 


In this section we discuss the geocoronal O l 1356 doublet, 
and the efforts which were undertaken to rule out emis¬ 
sion from this doublet in our G140L observation. At the 
redshift of M87, [Fexxi] has an expected wavelength of 
1359.9A, which is just a few resels away from geocoronal 
Ol A1358.51A. This is the weaker of the two lines in the 
Ol A1356 doublet; the stronger line has A = 1355.60 A. Ol 
A1356 dayglow emission is powered by collisional excitation 
from photoelectrons, while the nightglow is dominated by ra¬ 
diative recombination; the time-averaged dayglow emission 
looking towards nadir at 600 km (similar the orbital eleva¬ 
tion of HST, which is around 550 km) is more than an order 
of magnitude brighter than the nightglow ( Meier|1991 1. We 
use the timefilter routine provided in costools v. 1.0.6 to 
filter the data, keying on the altitude of the Sun in order 
to remove the data taken during local daytime (SUN_ALT 
> —10), which eliminates the dayglow emission. 

In order to be confident that we have accounted for 
all the geocoronal OlA1356, we briefly discuss the nature 
of the O iA1356 nightglow emission. Most of the O iA1356 
nightglow is produced at lower altitudes than the orbit of 
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